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ABSTRACT 



If dark matter in the form of compact objects comprises a large fraction 
of the mass of the universe, then gravitational lensing effects on gamma-ray 
bursts are expected. We utilize BATSE and Ulysses data to search for lenses 
of different mass ranges, which cause lensing in the milli, pico, and femto 
regimes. Null results are used to set weak limits on the cosmological abundance 
of compact objects in mass ranges from 1CT 16 to 1CT 9 M Q . A stronger limit 
is found for a much discussed Q = 0.15 universe dominated by black holes of 
masses ~ 1O 6 ' 5 M , which is ruled out at the ~ 90% confidence level. 



Subject headings: dark matter — gamma rays: bursts — gravitational lensing 
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1. Introduction 

The possibility that a cosmological abundance of dark compact objects (COs) in 
the universe could be revealed by gravitational lensing of more distant sources was first 
suggested by Press & Gunn (1973). The detection of lensed images from any source provides 
an estimation of the fraction of the universe composed of COs, and therefore the CO density. 
Non-detections offer a viable way to set constraints on the possible values of the CO mass, 
M co , and density, VL co . Depending on the angular separation of the lensed images, the 
gravitational lensing effect is classified into milli-, micro-, pico-, or femto-lensing. 

There is now overwhelming evidence that the majority of gamma-ray burst (GRB) 
sources lie at cosmological distances flMetzger et al. 1997] , [Kulkarni et al. 1998] , |Djorgovsk"i] 



st al. 19"9"B ). Their short duration and the transparency of the Universe in gamma rays 
make GRBs ideal probes of dark matter in the form of COs over the mass range 10~ 16 to 
1O 12 M . The probability of lensing a GRB depends on CO density (e.g. [Nemirofi 1989| , ), 
but it is estimated to be < 0.1% (piaes fc Webster 199| ). 



Multiple images of the same burst cannot be angularly resolved by present gamma-ray 
detectors, but they are easily temporally resolved (|Paczyhski 1987| ). Thus, rather than 
spatially identifying multiple images of lensed quasars, the lensed GRBs must be found 
by searching such events with identical temporal and spectral signatures. Marani et al. 
(1998) have compared the lightcurves of 1235 BATSE bursts, searching for gravitationally 
lensed GRBs caused by foreground galaxies, with negative results. The expected time delay 
between two detectable images can be written as (e.g. [Krauss fc Small 199 1|) 



where / is the brightness ratio of the two images also known as the dynamic range, is 
the lens redshift, R s ~ 3 km (M co /M Q ) is the Schwarzchild radius, and c is the speed of 
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light. Mao (1992) found that for a point mass lens, the median value of the time delay is 
At ~ 50 s x(M co /lO 6 M ). 

In this paper, we present results of searches for different types of gravitational lensing 
of GRBs, looking for any indication of the presence of dark matter in the form of compact 
objects. In all the monoluminous and non-evolving GRB source population is 

assumed. In §2, we describe the lensing types, search algorithms performed using BATSE 



data flMeegan et al. 19961) , and present results. We summarize and discuss these results in 
S3. 



2. Searching for Dark Compact Objects 
2.1. Millilensing 

In the case of millilensing, the COs have masses between 10 6 and 10 9 M & . The time 
delay between the images may vary from few seconds to ~ 10 4 seconds and their average 
angular separation is ~ 10~ 3 arcseconds. Note that M co = 1O 6,5 M is the Jean mass 
at recombination in many early universe scenarios and, consequently, a probable mass 
scale for dark matter. Gnedin & Ostriker (1992) proposed an open universe model with 
Qm — ^co — 0.15, filled with black holes of mass ~ 1O 6,5 M , which would dominate the 
dark matter galactic halos. 

Nemiroff et al. (1993) searched the first 44 BATSE bursts for lensed images that 
came shortly after the main event by autocorrelating each burst, and no lens candidates 
were found. This null result was interpreted by using the detection volume formalism of 
Nemiroff (1989), and by modeling the brightness distribution to estimate GRB redshifts 
( |Wickramasinghe et al. 1993]) . A standard (1,0) universe composed to closure density with 



objects of mass 10 6,5 - 10 8,1 M Q was marginally ruled out ([Nemiroff et al. 1993|) . 
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The data from Table I of Nemiroff et al. (1993) and their null results were used to set 
constraints on the CO mass and density for the Gnedin & Ostriker universe. Mao (1993) 
estimated that the lensing probability for this universe is ~ 5%. Following Nemiroff et al. 
(1993), for each burst we compute a detection volume between the observer and a burst at 
Zs, within which the lens with M co must fall to create a detectable image. Then, all 44 
volumes are added and the expected number of detectable lenses in our sample is obtained. 
For a given (Q co , M co ) pair, if the expected number is significantly greater than the null 
result obtained, we rule out such CO parameters. For a more detailed explanation of the 
method, see Nemiroff et al. (1993) and Marani (1998b). 

Figure |I] summarizes the results we obtain for the Gnedin & Ostriker universe. The 
figure shows the number of images expected for the 44 GRBs and for different (ft C o, M co ) 
parameters. The expected number of images for the Gnedin & Ostriker universe is ~ 2.5 - 3, 
in good agreement with Mao's estimations (1993): 2.2 lenses from 44 bursts were expected. 
Since no images were found, the universe with 

M co ~ 1O 6 - 5 M , Q M = Q CO = 0.15 (2) 

can be excluded at a confidence level of ~ 90.0%. 

2.2. Picolensing 

Picolensing refers to lensing mass range of ~ 10~ 12 -10~ 7 A'f©. The expected angular 
separation is ~ 10~ 12 arcseconds and the average time delay is ~ 10~ 18 seconds. The main 
lens candidates are primordial black holes, molecular clouds, and small planets. Currently, 
few arguments favor the hypothesis that extremely low mass objects compose a significant 
fraction of the dark matter. Nevertheless, if such objects were to exist, GRB lensing is an 
efficient method for their detection ([Nemiroff fc Gould 1995|) . 
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For these low mass COs, the corresponding radius of the Einstein ring projected onto 



our solar system, R E , will be at most 1 AU, Re = y (2R s D OL D os ) / D LS . Two spacecraft 
separated by 1 AU or more might see the same GRB event with very different amplification 
since, if the GRB were to lie within the Re of the lens for one spacecraft, it would lie outside 
Re for the other spacecraft. The i spacecraft observes the GRB with total amplification 



Ai = ((3 2 + 2)/(/3y (3 2 + 4), where j3 = (3/(3%, with (3 the angular separation between source 
and lens, and (3e the angular radius of the Einstein ring. There are currently two GRB 
detectors on orbit separated by a distance of at least 1 AU: BATSE and Ulysses GRB 
( [Hurley et al. 1992| ). If BATSE detects a strong burst while Ulysses determines it is a weak 



burst, the difference in the amplitude may be attributable to a gravitational lensing event 
due to masses in the range ~ 10~ 12 - 1O~ 7 M . 

We have compared the peak counts of 60 bright bursts which were triggered events in 
BATSE and Ulysses data (|Marani et al. 1998a| , [Marani 1998b| ). The energy range of Ulysses 



GRB detector is ~ 20-150 keV. The Ulysses triggered data have temporal resolutions of 8- 
and 32-ms, but they were rebinned into 64-ms count bins. For BATSE data, we have used 
64-ms concatenated ASCII data summed over energy channels 1 and 2, corresponding to 
the energy range ~ 20-100 keV. 

Once the data were background-subtracted and rebinned, the peak counts and the 
BATSE-to- Ulysses amplitude ratios, dA = C PiBATSE /C PiUlysses , were calculated. No evident 
discrepancy in the amplitude ratios has been found. Figure [2] shows BATSE peak counts 
versus Ulysses peak counts. The error bars are obtained assuming Poisson counts. The 
median value for the ratio, after normalization, is dA = 1, and the 2-cr upper and lower 
confidence levels are 2 and 0.68, respectively. A real difference in amplitude should be 
observable in all three timescales. We also compared peak counts for 256-ms and 1024-ms 
time resolutions, although no significant excess was found. 
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These null results are used to impose limits on the density and mass of COs, again 
using the detection volume formalism (|Nemiroff 1989|) . In the case of two observers, where a 



positive detection requires an amplitude difference between them, the volume within which 
a point-like lens must fall to create a discernible magnification is much more complicated 
than in the millilensing case, and no analytical expression is available. To find the expected 
number of lenses, we implemented a Monte Carlo simulation ( |Marani 1998b| ): For each 
burst in the sample, we randomly draw several lens positions on the lens plane - fixed at 
zl- We then compute the amplifications each observer would see if the lens is located in 
that particular position, take the ratio of those amplifications, and check if that ratio is 
above or below the 2-er C.L. Finally, we place the lens plane at another Zl, and re-compute 
the steps described above. When the lens plane reaches the source plane at zs, we obtain 
the detection volume V<i e t and the expected number of detectable lenses for that particular 
burst, given by Ndet = n co Vd e t- Here, n co = 3H%(1 + zl) 3 Qco/ (4ti-R s c 2 ) is the CO comoving 
number density. We then add all volumes of all bursts and repeat the calculations for 
different CO masses and densities. 

Figure shows the results we obtained with this simulation, assuming an Einstein-de 
Sitter universe. Less than ~ 0.9 detectable images are expected in our sample. Thus, the 
mass and density ranges of 

M co ~ 10" 12 - 5 - 1(T 9 - O M , n co > 0.9, (3) 

may be ruled out only at ~ 1-a statistical significance. 



2.3. Femtolensing 

The femtolensing mass range is M co ~1O" 16 -1O~ 13 M . The expected image separation 
is ~ 10~ 15 arcseconds, and the average time delay between the images is ~ 10~ 24 seconds. 
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Snowballs, primordial black holes, and comets are some examples of these COs. 

These extremely low mass lenses would produce interference fringes in the GRB energy 
spectrum. These features are evenly separated in energy space and might be confused with 
cyclotron absorption lines flGould 1992| ). Spectral lines between 10 and 100 keV were seen 
in missions like KONUS QMazets et al. 1982|) , HE AO 1 QHueter 19871) , and Ginga ([Fenimore 



et al. 1988| , |Murakami et al. 1988|) , and indeed have been attributed to cyclotron lines due 



to magnetic fields B ~ 10 12 G. However, these interpretations are still controversial. Briggs 
(1996) pointed out that many spectral line analyses performed in these GRB data sets used 
simplistic continuum models and incomplete statistical techniques. 

Several authors have analyzed BATSE spectra for spectral line features like the 



absorption lines detected by Ginga at ~ 20 and ~ 40 keV (see, for example, Band et al 



1996, Band et al. 1997, Briggs et al. 199£). Briggs et al. (1998) have recently performed a 



computerized search of 118 bright BATSE bursts. They have found only 4 candidates with 
one possible absorption line, but none of them with evenly spaced, multiple lines. Note that 
Stanek et al. (1993) and Ulmer & Goodman (1995) have argued that the detectability of 
the femtolensing signature depends on the source size. The more extended the source, the 
less pronounced the spectral lines. Lines at higher energies are progressively weaker for 
finite angular extent, since the source appears larger at higher energies. However, the short 
duration of GRBs (the vast majority lasting from milliseconds to a few tens of seconds) and 
their cosmological distances suggest that GRBs can be considered compact sources. 

Assuming that GRB sources subtend sufficiently small angular sizes to produce 
multiple, detectable lines in the spectra, the non-detections of spectral lines in 118 BATSE 
bursts by Briggs et al. (1998) can be used to impose lower limits on the mass and density 
of COs. Nemiroff (1989) estimated that the probability of lensing is ~ Q co /8 if the average 
distance to the GRBs is z GRB ~ 1, or ~ fi co /3 if zqrb ~ 2. Therefore, we can rule out, at 
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95 % confidence levels, universes composed of COs with M co ~ 10 16 — 10 13 M© and 

ttco > 0.2 , if z GRB ~ 1; or Vt co > 0.1 , if z GRB ~ 2. (4) 
The above estimations assume an Einstein-de Sitter cosmology and point mass lenses. 

3. Summary and Discussion 

The non-detections of lensing images in the current GRB data have been used, for 
the first time, to set conservative limits on the composition of dark compact object with 
masses between 10~ 16 and 10 -7 M & . Different searches in radio, optical, and gamma ray 
data have been conducted in the past to search for gravitationally lens-induced images 
caused by masses between 10~ 6 and 10 12 M Q ( |Carr 1994]) , but no search previous to the 
present one was carried out for the pico- and femto-lensing cases. The derived constraints 
on the density and mass parameters in different mass ranges are shown in Figure |j. The 
shaded areas show the excluded values of Q C o and M co . Light gray shadings represent 
very weak limits, whereas dark gray shadings are the imposed limits at ~ 2-cr confidence 
level. Three additional regions are shown in Figure f|. The cross-hatched area corresponds 
roughly to the range of values required to explain the MACHO and EROS microlensing 
results ( |Carr 1994]) . In the millilensing region, a universe proposed by Carr & Rees (1984) 



with (Qftj,Q\) = (1,0) and Q co ~ 0.15, M co ~ 10 6 M© is shown as a hatched, very small 
area. The constraints previously found by Nemiroff et al. (1993) are also illustrated. 
Using gravitational lensing data from the millilensing search done by these authors and 
the detection volume formalism ( [Nemiroff 1989| ), we were able to reject a popular universe 
proposed by Gnedin & Ostriker (1992) with M co ~ 1O 6 5 M and Q M = ^co = 0.15, at a 
confidence level of ~ 90%. 

No search for microlensing events (At ~ 5-10ms) in BATSE archival data was 
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conducted in the present work. However, Nemiroff et al. (1998) have argued that such short 
events might not be detected by BATSE due to its threshold setting for 64-ms timescale. 

The millilensing and picolensing limits are strong functions of the GRB redshift 
distribution, which is derived by modeling the observed brightness distribution. This model 
assumes no rate density evolution, standard candle, and a simple power-law spectrum, 
dN/dE oc E~ x . Several authors have recently suggested that the GRB rate density is 
proportional to the star formation rate in the universe QTotani 1997| , [Wijers et al. 199Sj) and 



that the intrinsic peak luminosities are described by a power-law distribution ( [Kommers 



|et al. 199~8"| |Krumholz et al. 1998| , [Mao fc Mo 1998| ). In both cases, larger redshifts are 



expected. Therefore, the probability of lensing by COs would increase and the limits 
imposed on fl co and M co would be stronger. 
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Fig. I. — Expected number of detectable lensing images for the given (Q co , M co ) values, in 
the (f2 M ,f2 A ) = (0.15,0) universe proposed by Gnedin & Ostriker (1992). The first available 
44 BATSE bursts are incorporated in this plot. 

Fig. 2. — Ulysses GRB peak counts versus BATSE peak counts for 60 bursts which were 
triggered events in both detectors. The solid line represents the median value for the 64-ms 
peak count ratio, whereas the dashed lines enclose the 2-a confidence level region centered 
in the median value. 

Fig. 3. — Expected number of detectable picolensing images for the given (Vt CO) M co ) values, 
incorporating 60 bursts that simultaneously triggered Ulysses GRB and BATSE. 

Fig. 4. — Constraints on the density and mass of individual compact objects, derived from 
the null results of gravitational lensing searches in BATSE data. 



